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ABSTRACT 

We estimate the contribution of Galactic pulsars, both ordinary and millisecond pul- 
sars (MSPs), to the high-energy (>100 MeV) 7— ray background. We pay particular 
attention to the high-latitude part of the background that could be confused with 
an extragalactic component in existing analyses that subtract a Galactic cosmic-ray 
model. Our pulsar population models are calibrated to the results of large-scale ra- 
dio surveys and we employ a simple empirical 7— ray luminosity calibration to the 
spin-down rate that provides a good fit to existing data. We find that while ordinary 
pulsars are expected to contribute only a fraction ~ 10~^ of the high-latitude 7— ray 
intensity {Ix ~ 1 x 10~^ ph s""'^ cm~^ sr~^), MSPs could provide a much larger con- 
tribution and even potentially overproduce it, depending on the model parameters. 
We explore these dependences using a range of MSP models as a guide to how 7— ray 
measurements can usefully constrain the MSP population. Existing 7— ray background 
measurements and source counts already rule out several models. Finally, we show how 
fluctuations in the 7— ray sky can be used to distinguish between different sources of 
the background. 

Key words: stars: neutron - pulsars: general - diffuse radiation - gamma-rays: 
theory 



1 INTRODUCTION 

Identifying the nature of diffuse radiation backgrounds 
is a fundamental problem because it may both lead to 
the discovery of previously unknown astrophysical sources 
and constrain the properties of known ones. While we 
now know the ultraviolet background to be dominated 



20091. High-energy (> 100 MeV) 7— rays are of particular 



by quasars and star-forming galaxies (e.g., Bolton et al. 
Faucher-Giguere et al.]|2008a|b[ |2009[ ) and that the 



2005 



soft X-ray background consists chiefly of emission from 
active galactic nuclei (AGN) (much of which has now been 
resolved by the Chandra X-ray Observatory and XMM New- 
ton; e.g., [Brandt fc Hasingerl|2005| |Hickox fc Markevitch| 



2006 20071, the sources of the diffuse 7— ray background 



have yet to be convincingly established. With the Fermi 
Gamma Ray Space Telescop^^ {Fermi) now taking data, it 
is timely to revisit this question from a theoretical point 
of view and make predictions that can be directly tested 
within the next few years. At the same time, the AGILI^ 
telescope of the Italian Space Agency is making interesting 
discoveries in spite of its more limited capabilities in the 
7— ray range (e.g., [Halpern et al.||2008| [Pellizzoni et al!] 
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interest as they probe some of the most energetic objects 
in the universe, involving sources such as pulsars (e.g.. 



Mukherjee et al. 



[Thompson et aL 1999 1 and blazars (e.g. 
'1997') emitting most of their electromagnetic luminosity in 
this band. 

In the standard picture, the 7— ray background is sep- 
arated into a strong Galactic component and a nearly 
isotropic component believed to be of extragalactic ori- 
gin. The Galactic component is generally thought to be 
predominantly produced by cosmic-ray interactions with 



interstellar gas and radiation (e.g., Bertsch et al. 1993 



Hunter et al. 1997 Strong et al.[[2000 l, while the isotropic 



Stecker 



component is commonly attributed to blazars (e.g 
fc Salamon|1996 Sreekumar et al.|1998 l. Some calculations 
however indicate that blazars may only account for a 
fraction of the total isotropic component (e.g., Chiang 



fc Mu kherjee] [19981 [Mukherjee fc Chiang[ [19991 [Miicke' 
fc Pohl 2000) and all blazar estimates presently suffer 



from significant uncertainties. For instance, the blazars 
seen by the Energetic Gamma Ray Experiment Telescop^ 
(EGRET) on the Compton Gamma Ray Observatory were 
almost always in a flaring state, with the duty cycle for 



http://heasarc.gsfc.nasa.gov/docs/cgro/egret 
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flaring being uncertain, and the evolution of blazar 7— ray 
emission must be extrapolated beyond the existing data 
I Thompson 2008 1 . Fermi has to date resolved only about 



7% of the isotropic component into blazars (Abdo et al 



2009b \ 



The separation of the 7— ray background into Galac- 
tic and isotropic components is a difficult task that relies 
on an accurate model for the Galactic contribution (e.g., 
Sreekumar et al.| [19981 |Keshet et aI1[2004l [Strong et aT] 
2004). It is therefore important to investigate Galactic 
sources that may contribute significantly to the high- 
latitude background and be confused with the extragalactic 
background. The detection of six pulsars at high confidence 
by EGRET ( [Thompson et al.|1999] ), the recent discovery of 
a radio-quiet 7— ray pulsar at the center of the supernova 
remnant CTA 1 ( Abdo et al.[[2008 l, and the firm detection 



of at least 29 pulsars in the 7— rays by Fermi (Abdo et al 



2009c I raise the possibility that pulsars could contribute a 
significant fraction of the high-latitude background. The 
high space velocities of pulsars (e.g. 



Gunn & Ostriker 



19701 [Lyne & Lorimer| 


1994| [Hansen & Phinney 


119971 


Arzoumanian et al. 


2002 


Faucher-Giguere & Kaspi 


[2006| 



imply that their spatial distribution may be extended above 
and below the Galactic plane and not correlate well with 
other tracers used to separate the Galactic component. 

Another timely motivation for understanding the sources of 
the 7— ray background is provided by indirect searches for 
dark matter annihilation products, which may bo detectable 



in 7— raj 


s (e.g., 
al. 200 


Diemand et al 


[2007 Kuhlen et al., 2008( 


Baltz et 


8; [Springel et al. 2008| [Kamionkowski &[ 


Profumo 


|2008( |Dodelson et al. 


2009 


I. Among the possible 



astrophysical sources, 7— ray pulsars might be the most dif- 
ficult to disentangle from a dark matter annihilation signal 
owing to their similar spectra ( [Baltz et al.|2007[ ). Given the 
low intensity expected from dark matter annihilation, even 
a small contribution by pulsars to the 7— ray sky could be 
an important contaminant for these searches. 

The idea of pulsars as contributors to the 7— ray background 
has previously received some attention in the literature. Or- 
dinary pulsars have been suggested to provide a significant 
fraction of the low-latitude Galactic 7— ray background, 
with estimates ranging from 5-10% to comparable to the 



cosmic ray contribution ( 


ligdon & Lingenfelter 


1976 


Strong et al.|1977| Harding & Stecker|1981 [Harding 


1981a 


Bhattacharya & SrinivasanI 1991[ [Sturner & Dermer|[1996 


McLaughlin & Cordes 200C 


) 1 . It has also been proposed 



that millisecond pulsars (MSPs), owing to their rapid spin 
frequencies, may be bright in the 7— rays ( [Usovj [1983| 



Sturner fc DermCT][1994 Harding et al.|2005 I and therefore 
also contribute significantly to the 7— ray background 
([Bhattacharya fc Srinivasan[[T99T| [Bailes fc Kniffen|[T992[ 
Bhatia et al.[[1997[). 



The present work aims to improve previous crude es- 
timates of the pulsar contribution to the high-energy 7— ray 
background by using the pulsar population synthesis code 
developed by [Faucher-Giguere fc Kaspi] ( [2056| ) (FGK06), 
with the goals of making predictions for upcoming Fermi 
observations and providing a guide to use these to constrain 



the population and astrophysics of Galactic pulsars. The 
population synthesis code allows us to consider pulsar pop- 
ulations with realistic distributions of physical properties 
and calibrated against existing large-scale radio surveys. 
As we will argue, however, simple physical arguments 
and comparison with existing 7— ray data also allow 
relatively model-independent conclusions to be drawn. 
We pay particular attention to the pulsar contribution to 
the high-latitude background that may be confused with 
the extragalactic component or the dark matter signal. A 



similar approach was used by Strong ( 2007 1 to quantify the 



contribution of unresolved sources to the Galactic diffuse 
emission. 

We begin by describing the population synthesis mod- 
els used to calculate the contribution of pulsars to the 
7— ray background in i|2] Our results are compared with 
existing 7— ray observations in 33| and the implications 
for pulsar physics are discussed ^ In !|5] we show how 
fiuctuations in the 7— ray sky can be used to discriminate 
between different origins of the 7— ray background. We 
finally conclude in i|6] 



2 POPULATION SYNTHESIS MODELS 

We first estimate the expected pulsar contribution to the 
7— ray sky by combining pulsar population synthesis mod- 
els calibrated against radio surveys with a prescription mo- 
tivated by existing 7— ray observations for the 7— ray lumi- 
nosity of pulsars. 



2.1 Models 

The basis for our calculations of the 7— ray emission from 
pulsars is the population synthesis code developed by 
FGK06, who studied the birth properties and evolution 
of isolated radio pulsars. In that work, a Milky Way 
model was seeded with nascent pulsars which were phys- 
ically evolved in time. The synthetic pulsar population 
was then observed by mock realizations of the Parkes 



and Swinburne Multibeam Pulsar Surveys (Manchester 
[et al.[[200 r: Edwar ds et al.[[200l| ), including all the main 
selection effects, and the results compared with the ac- 
tual detections. This resulted in a model for the spatial, 
kinematic, rotational, magnetic, and radio properties of 
the pulsar population that reproduces the observations well. 

In this work, we use this model to calculate the 7— ray 
contribution of ordinary pulsars. As will be discussed below, 
we extend the code to also model a population of MSPs. 
The evolution of these recycled pulsars, whose spin-up to 
millisecond periods is thought to have occurred by angular 
momentum accretion from a binary companion, is much 
more complex (e.g., [Bhattacharya fc van den Heuvel|199T| ). 
We therefore do not attempt to model their evolution but 
rather prescribe probability distributions for the properties 
of their current Galactic population. This population is 
also not as well constrained observationally, both due to 
the relatively small number of known MSPs and to the 
difficulty of accurately modeling the selection effects that 
affect their detection. For example, the scattering of radio 
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waves by free electron density fluctuations (e.g., Rickett distribution (Cordes & Lazio 20021 given the sightline 



19901 makes it challenging to detect short pulsar periods 



through the Galaxy. The sky brightness temperature at 408 



toward the Galactic center and the orbital acceleration in 



binary systems acts to mask periodicity (e.g., Johnston & 
Kulkarni 19911. Another potentially important selection 



effect against the detection of MSPs is that obscuration 
of the pulsar signal by material blown off a companion 
by the pulsar wind may render it invisible, especially in 



the radio (Tavani 19911. The fastest-spinning known MSP 
J1748— 2446ad is in fact eclipsed in the radio for ~ 40% of 
its orbit ( [Hessels et"aL]|2006 K 



As will be argued (see [2.8 and however, simple 

physical arguments can be used to draw conclusions on 
the contribution of MSPs to the high-latitude 7— ray back- 
ground that are insensitive to the details of the population 
models (see also Strong 2007 for a similar approach in the 
context of the Galactic diffuse emission). 



2.2 Calibration Against Radio Surveys 

Most of the constraints on pulsar populations come from 
radio observations, with about 1800 objects now docu- 
mented in the ATNF Pulsar Database|fJ Since we will 
explore models beyond the ones studied by FGK06 for 
MSPs, we pause to summarize how we model the sensitivity 
of radio surveys to pulsars. The mock surveys will be used 
to calibrate our models based on the actual observations. 

The detectability of a pulsar depends on its intrinsic 
properties (luminosity, pulse period, and duty cycle), 
on its location (distance, dispersion measure, interstellar 
scattering, and brightness temperature of the background 
sky), and on the details of the observing system. The 
minimum flux theoretically detectable from a radio pulsar 
is estimated by the usual radiometer equation, 



'S'min — 



/^^sys'^(7rcc ~l~ ^sky) 



P-We 



(1) 



where 5bcam is a factor accounting for the reduction in 
sensitivity to pulsars located away from the center of the 
telescope beam, Tree is the receiver temperature on a cold 
sky, Tsky is the sky background temperature, G is the 
antenna gain, A'p is the number of polarizations summed, 
Af is the receiver bandwidth, tint is the integration time, P 
is the pulse period. We is the effective pulse width, a is the 
signal-to-noise detection threshold, and /3aya is a constant 
accounting for various system losses (e.g., Dewey et al. 
1984|. 



Further details regarding the 5min calculation are pro- 
vided in Appendix 1 of FGK06. Briefly, the effective pulse 
width Wa is a quadratic sum of the pulsar's intrinsic 
pulse width, the effective sampling time of the survey, the 
dispersion smearing across a receiver channel, and the pulse 
scattering time scale. The dispersion measure and pulse 
scattering time of each simulated pulsar are calculated using 
the NE2001 model for the Galactic free electron density 



http:/ /www. atnf.csiro.au/research/pulsar/psrcat/ i jManch-" 
es ter et aI.|2Q05| . 



MHz is taken from the Haslam et al. ( 1981 1 all-sky map and 
scaled to the desired freque ncy assuming a spe ctral index 
Qbg = 2.7 (rsky(f) oc i/~">'g; |La.wson et al.|l987| . The 5bcam 



flux degradation factor is calculated assuming a Gaussian 
power pattern, with each synthetic pulsar in the covered 
sky area lying at a random position within the half-power 
cone of a telescope beam. A synthetic pulsar is detected if 
and only if it lies in the sky area covered by a survey, the 
radio flux density of the pulsar at the central frequency 
of the survey at the Earth exceeds the survey threshold 
given by equation ([TJ, and the pulsar is actually beamed 
toward the observer in the radio. The latter criterion is 
implemented by assigning a probability /b'"^'" (which in 
general can depend on the period of the object) to each 
pulsar that its beam crosses our line of sight to it. 

The radiometer equation |l| gives the minimum flux 
for the detection of a pulsed signal and is valid for solitary 
pulsars. While less than 1% of ordinary pulsars are observed 
to have orbiting companions, about 80% of MSPs are found 
in binary systems (e.g., Lorimer 20051. The periods of 



pulsars in binary systems are modulated by the Doppler 
effect induced by orbital acceleration. This acceleration 
tends to mask the periodicity of the signal by spreading it 
into multiple Fourier bins (e.g., Johnston fc Kulkarni|199T |. 
The resulting degradation in sensitivity is a function of the 
orbital parameters of the binary system and is in general 
most important when the orbital period is less than the 
integration timej^The most disfavored pulsars are therefore 
MSPs with close and/or massive companions. As a result, it 
is difficult to quantify the number and characteristics of the 
binary pulsars that are missed in existing surveys. Most, if 
not all, MSP population studies to date have neglected these 
effects. In this work, we simply note that the unmodelled 
degradation in sensitivity to radio MSPs in binary systems 
can only lead us to underestimate the number of MSPs in 
the Galaxy (by optimistically estimating their detection 
probability) and thus to underestimate their cumulative 
7— ray emission in particular models. 



2.3 Gamma-Ray Luminosity Prescription 

In order to calculate the pulsar contribution to the 7— ray 
background, we must assign a 7— ray luminosity 
to each pulsar in our radio-calibrated populations. No 
definitive prescription for doing so exists at present and 
Fermi is poised to provide unprecedented observational 
constraints on theoretical models. Our approach is to adopt 
an empirically-calibrated prescription, which we will use 
to make predictions that can be tested with upcoming 
observations. 

Unlike in the radio, where the emitted power is a negligible 
fraction of the total emission, a significant fraction of the 
rate of loss of rotational kinetic energy, E, of pulsars is 



^ For integration times that arc much less than the orbital period, 
the Doppler effect is nearly constant during the observation and 
the periodicity is approximately preserved. 
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observed in the form of 7— rays. In fact, the maximum 
observed energy output is in the 7— ray band for all the 
pulsars detected by EGRET. For these sources, the ratio of 
the electromagnetic power above a photon energy of 1 eV to 
the total spin-down power ranges between about 0.1% and 
20%, Eissuming a 7— ray beam covering 1 steradian (Rudak 



fc Dyks||1998[ [Thompson et aLH iggQ*). For a fiat spectrum, 



with equal power per energy decade, half of this power is 
in the form of > 100 MeV 7— rays. The flux measured by 
Fermi for the pulsar at the center of the CTA 1 supernova 
remnant also indicates a 7— ray efficiency of about 1% for 
the same beaming solid angle, but perhaps as high as 10% 
if the emission arises in the outer magnetosphere, for which 
models predict beam sizes ^ 1 sr ( Abdo et al.||2008 l. 



The few pulsars studied in the 7— rays by EGRET 
suggest a simple relation 



E oc py^p-y^^ (2) 

where E = 4tv^IP/P^ is the rate of loss of rotational kinetic 
energy of the neutron star, with I (conventionally assumed 
equal to 10*^ g cm^) denoting its moment of inertia. This 
relation was first noticed and theoretically investigated by 
Arons (19961 (see also Thompson et al. 19991. No other 



simple parameters appear well correlated with the observed 
properties of these pulsars ( [Goldoni fc Musso| |1996| . In 
particular, there appears to be more variance in the fraction 
of E that is converted into 7— rays. 

The empirical relation in equation Q can also be motivated 
theoretically by the facts th at the total voltage drop in the 
magnetosphere AV oc 



E ([Ruderman fc Sutherland] 1975 

Arons|1996 1 or that the Goldreich & Julian (19691 current. 



defined as the rate of outflowing relat ivistic particles over 
the polar cap, A'^gj oa 



E as well (Harding 



1981b I. In 



outer gap models of high-energy emission, this relation 
can also be approximately derived as the product of the 
volume of the outer gap, the particle number density 



assumed to be close to the Goldreich & Julian ( 1969 1 value. 



and the power radiated by curvature radiation by a single 
particle (e.g., Romani|1996 Crusius-Waetzel fc Lesch|2002 |. 



To avoid spectral uncertainties in the conversion from 
energy units to photon units, we express luminosities 
directly in units of ph s^^ above 100 MeV, and denote 
them L^^. Since the 7— ray fluxes reported by Fermi are 
already expressed in ph s~^ cm^, our empirical calibration 
( ^2.4[ ) is free of any assumption with regards to the slope 
or high-energy cutoff of the 7— ray pulsar spectrum. More 
precisely, we define the 7— ray photon luminosity of a pulsar 
above 100 MeV to be given by 



= Kmm{CP'^'^p-'"\ (3) 

where C is an empirically-determined proportionality con- 
stant, Z™^" is the maximum fraction of the rotational power 
loss converted into 7— rays, and K accounts for the conver- 
sion from erg s~^ to ph s~^. 



2.4 Observed Gamma-Ray Luminosities 

The parameters of our 7— ray luminosity prescription are 
set based on the measured properties of the known 7— ray 
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Figure 1. Estimated 7— ray photon luminosity at ener- 
gies >100 MeV vs. rate of loss of rotational kinetic en- 
ergy for the firm pulsar detections in the first three months 
of Fermi data and in the first year of AGILE data for 
which rotational parameters and distances are tabulated in 
the ATNF database. The crosses show the EGRET pul- 
sars (from left to right: B1055-52, B1951+32, B1706-44, 
J0633-I-1746, B0833-45, B0531+21). The triangles indicate 
MSPs (J0030-I-0451, J2124-3358, J0613-0200, B1821-24*). The 
filled circles show other ordinary pulsars detected by Fermi 
or AGILE (J0631+1036, J1509-5850, J1028-5819, B1046-58, 
J2021-I-3651, B1509-58*, J2229-I-6114*). The asterisks indicate 
AGILE discoveries. The point for B1046— 58 was translated to 
the left by AE = — 5 X 10~^^ erg s~^ to distinguish it from MSP 
B1821-24. The soUd line shows the best-fit LP*" oc E law and the 

E law. The MSP points were ex- 
cluded for these fits. The dotted line indicates L-y = 0.05-E under 
the assumption of a spectrum with equal power per decade energy 
and a high-energy cutoff -Bmax = 3 GeV. The 7— ray luminosity 
law adopted in our calculations is shown by the broken thick grey 
line. The calculated 7— ray luminosities assume isotropic emis- 
sion. 



pulsars. Specifically, we consider the subset of the 29 
firm pulsar detections in the first three months of Fermi 
data ( Abdo et al.||200"9c l whose rotational parameters and 
distances are tabulated in the ATNF database and the 
three firm first-year 7— ray discoveries by AGILE (two of 
these are not part of the firm Fermi detections; |Pellizzonil 
20091. This excludes 14 objects newly discovered by 



dashed line the best-fit L^^ 



et al. 



Fermi and for which these parameters (usually derived 
from radio counterparts, which may not exist in these 
cases) are not yet available. This however includes all the 
EGRET pulsars and three 7— ray MSPs. Figure [T] shows 
the 7— ray photon luminosity for energies above 100 MeV 
for these pulsars versus their inferred spin-down luminosity. 
The 7— ray photon luminosity was in each case estimated as 
= AiTcfFP'', where F^"" = F23 + F35, and F23 and F35 
are the photon fluxes for energies 100 MeV— 1 GeV and 1 
GeV- 100 GeV for each pulsar given by [Abdo et aL] ( |2009c| ). 
For the AGILE detections with no Fermi counterpart in 
Abdo et ak] (|2009c|l , the photon fluxes at >100 MeV were 
taken from |Pellizzoni et al.| ( |2009"| ). The spin-down luminosi- 
ties E and best distance estimates d were taken from the 
ATNF database. One exception is for PSR J2021-I-3651, for 
which we assumed a distance 12 kpc quoted by |Hessels et al.| 
( |2004| ) and consistent with the NE2001 free electron density 
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model. The error bars account for a factor of 2 uncertainty 
in d, roughly what can be conservatively expected for 
distances determined from the dispersion measure. Variance 
in the viewing geometry may introduce additional scatter. 
The uncertainty in the 7— ray flux is generally subdominant. 

Figure [l] also shows best-fit L'^^ cx E and oc 
lines. For these fits, the MSPs were excluded as a test 
of how well the fits to the better-studied ordinary pulsar 
data predict the MSP 7— ray luminosities. The prescrip- 
tion L^^ oc \/S suggested by the EGRET pulsars holds 
reasonably well for this extended data set if we take into 
account the intrinsic scatter that may exist in the true 
7— ray luminosity law and the uncertainties in the distances 
and beaming angles of the pulsars. The extrapolation to the 



firm pulsar detections in the first three months of Fermi 
data are new discoveries with no previously known radio 
counterparts. Since it is in principle easier to detect and 
confirm 7— ray pulsars associated with known sources, this 
suggests that the 7— ray beaming fraction /J is in general 
much larger the radio beaming fraction /b'''*'°, and likely 
of order unity. Theoretical models in which the 7— ray 
emission originates from the outer or slot gaps of pulsar 
magnetospheres in fact predict large 7— ray beams of size 
1 sr (e.g., Romani & Yadigaroglu 1995 Hirotani 2008 



Harding et al. 20081. Comparatively large 7— ray beaming 



fractions would also be consistent with many of the EGRET 
unidentified sources being pulsars in spite of the mostly 
unsuccessful searches for radio pulsars within their error 



boxes (e.g., Kramer et al. 2003 Champion et al. 2005 



MSPs is excellent, in agreement with |Abdo et al.| (|2009a[), [Crawford et al.|2006| ). Note that in the limit of large 7-ray 



who find that the basic emission mechanism seems to be 
the same for MSPs and young pulsars. Again, because the 
empirical luminosities are already expressed in ph s~^, the 
resulting best-fit relations for L?i^ = L?f'^{E) do not depend 
on the spectral parameters. While this has no impact at 
all on our other calculations, we may assume that pulsars 
have spectra containing approximately equal power per 
energy decade up to a high-energy cutoff -Emax ~ 3 GeV to 
convert thephoton luminosities to cgs units. The dotted 
line in Fig. Ill indicates Lj = 0.Q5-E for this fiducial choice 
and shows that physically limiting the 7— ray luminosity 
at 0.05E (Z™'^" = 0.05) provides an excellent match to the 
MSPs. The resulting 7— ray luminosity law, used in the rest 
of this work, is shown by the broken thick grey line. 

After our analysis was completed, the first Fermi 7— ray 
pulsar catalog, containing more pulsars than we have 
analyzed, was released (Abdo et al. 2009c I. We have 
verified that our best-fit LJJ" = L?j'^{E) relation is broadly 
consistent with the expanded data set. Moreover, as our 
main results are only weakly dependent on the details of 
the population synthesis assumptions by virtue of their 
calibration to existing source counts (see [2.8 and they 
are unaffected by the new data. 



2.5 Gamma- Ray Beaming Fraction 

The 7— ray luminosities above were derived assuming a 
7— ray beam size of 47r sr. This assumption has no effect on 
the 7— ray flux observed at the Earth for any given pulsar; 
while the inferred 7— ray luminosity of the pulsar scales 
with the beam size, so does the factor by which the flux is 
diluted as it propagates toward us. However, larger beam 
sizes imply that more pulsars have their beams intersecting 
our sightline to them and are actually observable from the 
Earth. This effect therefore enhances the predicted 7— ray 
background for particular models. Because the shape of 
the flux distribution is determined predominantly by the 
geometry of the spatial distribution (see models which 
predict the same number of detections above a given flux 
threshold also however predict approximately the same 
integrated 7— ray background. 



beams, the details of the beaming geometry have no impact 
on predictions of the integrated 7— ray background. 

The beaming fractions of ordinary and millisecond 
pulsars in the radio, which are important in normalizing 
their populations to radio survey detections, will be treated 
separately in §2.6|and §2.7[ 



2.6 Ordinary Pulsars 

For the ordinary pulsars, the population parameters 
and number normalization were determined in detail by 
FGK06 through a comparison with the 1.4 GHz Parkes 
and Swinburne Multibeam Pulsar Surveys ( jManchesterj 



et al. 2001 Edwards et al. 20011. These two surveys, 



which together detected over 1000 of the 1800 pulsars in 
the ATNF database and which were performed with the 
same observing system, provide a comparison sample with 
a homogeneous selection function. The ordinary pulsar 
population parameters are therefore taken as fixed. 

Briefly, the synthetic population of ordinary pulsars is 
an evolved snapshot of pulsars that are born in the Galac- 
tic spiral arms with a Galactocentric radial distribution 
peaking at radius ~ 3 kpc and with a mean birth velocity 
of 380 km s~^. The birth spin period distribution is normal, 
with (Po) = 300 ms and apg = 150 ms, truncated so that 
Po > 0. The magnetic field distribution is lognormal, with 
(logB) = 12.65 and (Jiogs = 0.55. Throughout, magnetic 
fields are expressed in G and log refers to the base- 10 
logarithm. The pulsars are evolved in time by solving for 
their orbits in the Galactic potential and by assuming 
magnetic dipole breaking with constant magnetic fields. 
The radio luminosities of pulsars are also assumed to be 
proportional to as was found to produce a good fit to 

the P — P diagram, but with a random dithering factor. 
The period-dependent radio beaming fraction is the one 



obtained by Tauris & Manchester ( 1998 1 in their analysis 



of polarization data: 

/radio, ord / tj\ 
b 



0.09[log(P/s) - 1]^ + 0.03. 



(4) 



2.7 Millisecond Pulsars 



The beam sizes of pulsars determine the fraction /b 
of all pulsars whose beams cross our sightline. Half of the 



Our population synthesis treatment for MSPs differs from 
our treatment of ordinary pulsars in that we do not attempt 
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to model their time evolution. Rather, we populate a 
present-day synthetic Galaxy with MSPs whose properties 
are drawn from prescribed distributions. Since the con- 
straints on this population are poorer, due to their more 
complex evolutionary channels and more severe selection 
effects, we will explore a range of models. This approach 
will also serve as a guide of how 7— ray observations can 
provide new constraints on MSPs. In this work, we do not 
attempt to constrain the details of the MSP population 
with radio data as was done for ordinary pulsars in FGK06 



(for work in this direction, see Cordes & Chernoff 1997 



|Lyne et aL]|1998| [Story et ar]|2007| |. Rather, the present 



simulations are simply meant to illustrate the potential 
contribution of MSPs to the 7— ray background based on 
fiducial population parameters in the literature. As will be 
shown, rather model-independent conclusions can be drawn 
on the contribution of MSPs to the 7— ray background 
and the most accurate estimates are likely based not on 
constraining the details of the population parameters from 
radio data, but by calibrating to existing 7— ray results 
(^. The parameters we explore are tabulated in Table [5] 
and the model framework is described below. 

As a starting point, our models are motivated by the 
MSP population synthesis study of [Cordes fc Chernoff| 
(19971. Another useful reference is the population analysis 



of the 400 MHz Parkes Southern Pulsar Survey ( Lyne et al 



1998). Story et al. (20071 also recently performed a MSP 



population synthesis study. The important MSP parameters 
for our purpose are their spatial distribution and their 
distributions of period, magnetic field, and radio luminosity. 

Before proceeding, we pause to comment on the bi- 
nary nature of most MSPs and their subset associated with 
globular clusters (GCs). As discussed in [2.2 about 80% 



of MSPs have a binary companion. We do not explicitly 
model binary systems in this work. Rather, we assume the 
distributions we prescribe to apply to the MSP population 
as a whole, including those in binary systems. The prin- 
cipal effect that is missed in this approach is the reduced 
sensitivity of surveys to pulsars that are either accelerated 
or eclipsed in binaries. As noted in §2.2| the unmodelled 
degradation of the sensitivity to radio MSPs in binary 
systems leads us to optimistically estimate the sensitivity 
to MSPs in our simulations and therefore to underestimate 
the number of MSPs in the Galaxy and their contribution 
to the 7— ray sky. While individual GCs can harbor large 
concentrations of MSPs, the total number of GC MSPs 
in the Galaxy is much smaller than the total number of 
field MSPs. GCs will therefore appear as bright 7— ray 
point sources, but will contribute negligibly to the 7— ray 
background. We therefore do not model GC MSPs here. 

FGK06 found the best match to the observed Galactic lon- 
gitude distribution of ordinary pulsars for a Galactocentric 
birth radial distribution like peaking at « 3 kpc from the 
Galactic center (like the one advocated by |Yusifov fc KiigiikI 
2004 1 . Although this is consistent with the distribution 



of the massive stellar progenitors of pulsars, MSPs have 
characteristic spin-down times up to 10 Gyr or more (e.g.. 



Camilo et al. 1994 1 and so can be generally expected to 



Their current positions should therefore be largely uncorre- 
lated with their birth locations. In particular, the annular 
overdensity and spiral arm structure should be almost 
completely smeared out. Thus, for MSPs, we adopt a model 
in which the surface density projected on the Galactic plane 
is described by 



p(r) cx exp (— r /2cr^); < r < 00, 



(5) 



where r is the distance from the Galactic center. Based on 
the birth distribution of ordinary pulsars, a value (Jr ~ 5 
kpc is expected. 

The height distribution above and below the Galactic 
plane is taken to follow a simple exponential distribution. 



N{z) oc exp { — \z\/{\z\)), 



(6) 
plane 



parameterized by the mean distance from the 
{\z\). For this model, [Cordes fc Chernoff| ( pff]) found 
(l^l) — 0.5 kpc. Using dynamical simulations, Lyne et al. 



(19981 found that the underlying \z\ distribution can 



be quite broad, with a median of 570 pc but a mean 
of 1.1 kpc, if the velocity distribution is a Maxwellian 
with a mean 3-D velocity of 130 km s"'^. With similar 
simulations, but with a mean 3-D velocity of 110 km 
s~^ and a different model for the Galactic gravitational 
potential. Story et al. ('2007) found a \z\ scale height of 
0.5 kpc, but also with indications of an extended large \z\ 
tail. We will explore different values from (z) — 0.5 to 2 kpc. 



Lorimer et al. ( 1996 1 studied the underlying spin pe- 



riod distribution of MSPs and highlighted the difference 
with respect to the observed distribution stemming from 
the selection effects against detecting the shorter period 
pulsars. In their likelihood analysis of these selection effects. 



Cordes & Chernoff (19971 assumed a power-law functional 



form for the underlying spin period distribution. 



N{P) oc^ D-P.2-1 



pmpi- 
P" 



if -Pmin 4, X < 
if P ^ Pbroak 



(7) 



have gone through many orbits in the Galactic potential. 



and found a best fit for mp,i = mp^2 = ^1, with Pmin > 0.65 
ms at 99% confidence and Pmin > 1 ms at 95% confidence. In 
their analysis a single slope parameter mp = mp,i = mp,2 
was assumed; we generalize the functional form here to allow 
for a fiattening at short periods. The parameters mp^i and 
mp,2 correspond to slopes in a log A^-log P diagram. The 
fastest spinning pulsar known at present, J1748— 2446ad, has 
a period of 1.4 ms ( Hessels et al.|2006 l, constraining Pmin to 
be approximately below this value. Similarly to the ordinary 
pulsars, we model the magnetic field distribution of MSPs 
using a lognormal distribution, 

iV(logB) (X exph(logB - (logP))V2afogB], (8) 

with mean logarithm (log B) and logarithmic standard devi- 
ation CTiogB- The period derivative, P, is then related to the 
period and magnetic field assuming magnetic dipole braking 
using the conventional formula B = 3.2 x 10^^(PP/s)^/^ G 
(e.g., [Manchester fc Taylor||1977| ). 

The radio luminosities of the MSPs, conventionally 
defined as Lr = Srd^ (where St is the fiux density at the 



radio frequency of interest), have been inferred by Lyne 
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et al. ( 1998 1 to be well described (after correcting for 



selection effects) by a log A'^-log Lr distribution of the form 



Tvr/ r \ I -^r -t^r.min ^ I^r \ I^r, break /n\ 

N{L,) oc <; ™^ ,_i , _ . (9) 

i^r n _ijr ^ ^r, break 



with m_L = rriL,! = mL,2 = —1, in agreement with ordinary 



pulsars. Cordes & Chernoff ( 1997 1 used the same functional 



form and found the same slope in their likelihood analysis. 
We again have generalized the distribution to allow for a 
change of slope at Lr, break- Contrary to the flattening that 
is observed for ordinary pulsars below = 20 mjy kpc^, 
Lyne et al. ([l998) found a slope ttil = — 1 down to the 
lowest luminosities sampled ~ 0.3 mJy kpc^ at 400 MHz, 
suggesting that the number of low-luminosity radio MSPs 
could exceed that of ordinary pulsars. 

For consistency with the ordinary pulsar analysis of 
FGK06, the radio luminosities of the MSPs in our sim- 
ulations are set at 1.4 GHz. To convert to luminosities 
at arbitrary frequencies, each MSP is assigned a spectral 
index (S,^ oc from a normal distribution with a mean 

(a) = 1.8 and a standard deviation Ua = 0.55 (Kramer 
et al.|[l998 l. Like ordinary pulsars, the radio emission from 
MSPs is assumed to cease when they cross the "death line" 
in the P — P diagram (eq. 22 in FGK06). Because of the 
short periods of MSPs and since the fastest-spinning ones 
are assumed to be the strongest 7— ray emitters, the death 
line plays a minor role in regulating the cumulative 7— ray 
emission from MSPs in our models. The beaming fractions 
of MSPs tend to be lower than predicted by fits to the 
longer-period pulsar population, as in equation Q (e.g., 
[Kramer et al.|[T998| ), which would imply beaming fractions 
near unity. We fiducially take the beaming fraction of MSPs 
to be j^'"i'°''^SP _ g though our results can be easily 
scaled to any other constant value (eq. |13[ ). 

Of all the pulsar surveys (excluding deep searches to- 
ward globular clusters), the MSP detections are dominated 
by three. They are the Parkes and Swinburne Multibeam 
Pulsar Surveys that were also used to constrain the popu- 
lation of ordinary pulsars. In addition, the Parkes Southern 
Pulsar Survey which covered all southern hemisphere 



declinations at 400 MHz ( Manchester et al. 1996 1 took 



advantage of the steep spectra of pulsars in the radio and 
was sensitive to low-luminosity pulsars and MSPs. These 
surveys detected 20, 21, and 19 MSPs, respectively. Our 
MSP population simulations use 2 x 10^ Monte Carlo 
pulsars each. Quantities that depend on the total number 
of pulsars in the Galaxy, such as the 7— ray background 
they produce, are calculated by normalizing the number 
of pulsars to radio detections. Let A'^^'^'*' and A'^l"" be the 
numbers of real and simulated detections for a survey 
labeled x. We use a normalization weighted by the numbers 
of detections in the real surveys: 



(10) 



where x runs over the Parkes and Swinburne Multibeam 
surveys and the Parkes Southern survey. 



2.8 Ordinary vs. Millisecond Pulsars 

In this section, we outline astrophysical considerations 
that will be helpful in understanding which population, 
the ordinary pulsars or the MSPs, may be expected to 
have the greatest contribution to the high-latitude 7— ray 
background. Some of these considerations will also allow 
us to draw robust, model-independent conclusions with 
respect to the 7— ray background contribution of MSPs, for 
which the population parameters are most uncertain. 

Let us first focus on the 7— ray luminosity of individ- 
ual sources. Assuming that the 7— ray luminosities of both 
populations follow the simple relation in equation ([3|, we 
can express the ratio of the 7— ray luminosity of a MSP to 
that of an ordinary pulsar. 



r ord 

^1 



1/2 



-3/2 



(11) 



where the subscripts MSP and ord refer to values char- 
acteristic of MSPs and ordinary pulsars, respectively. In 
this expression, we have neglected the fact that the 7— ray 
luminosity is physically limited to a fraction of the spin- 
down luminosity of the pulsar. Taking values Pord = 0.5 



s, Pm 



3 mS, Pord 



10" 



and Py 



10" 



s , roughly at the center of the observed populations in 
the P - P diagram, we find L^^^ / L°^"^ ~ 20. As a result 
of their short periods, a "typical" MSP may therefore be 
brighter in the 7— rays than a "typical" ordinary pulsar. 

Since L^y oc \/^ in this model, the conclusion remains 
qualitatively true even if the 7— ray luminosity is limited at 
a fraction of E. 

Note that the sample of observed ordinary pulsars 
may be biased toward short-period pulsars if young objects 
tend to be brighter in the radio (e.g., FGK06). On the 
other hand, the most rapidly rotating MSPs are strongly 
disfavored in existing surveys. It is therefore likely that 
the known ordinary pulsars are biased toward the 7— ray 
bright end of their population, whereas the known MSPs 
may be biased toward their 7— ray faint end. As a result, 
the truly typical MSP may be brighter in the 7— rays by a 
substantially larger factor than the truly typical ordinary 
pulsar than estimated above. 

The above argument is perfectly consistent with the 
few observed ordinary pulsars to appear more luminous in 
the 7— rays than the few observed MSPs, as suggested by 
Figure [1] In fact, the brightest members of these respective 
populations are not representative of their typical members. 
Any meaningful comparison must take into account the 
very different luminosity distributions of these populations. 
We have explicitly verified using our population models that 
the 7— ray luminosities of the ordinary pulsars predicted to 
have been detected in the early Fermi data are generally 
higher than the same quantity for the detected MSPs, in 
agreement with the data shown in Figure [l] At the same 
time, the average member of the underlying MSP popu- 
lation is predicted to be more luminous than the average 
member of the underlying ordinary pulsar population, and 
MSPs are predicted to provide a higher contribution to the 
7— ray background. This finding strongly cautions against 
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drawing conclusions from the observed "tip of the iceberg" 
of the 7— ray pulsars. 

How many Galactic MSPs are there compared to or- 
dinary pulsars? 

If MSPs are recycled pulsars, then averaged over a 
long time scale the birthrate of MSPs must necessarily 
be less than the birthrate of ordinary pulsars, which may 
viewed as the rate of "pulsar creation." Based on this 
consideration only, the MSP birthrate can be expected to 
be substantially lower than that of ordinary pulsars since 
only pulsars in binary systems with specific properties (like 
low-mass X-ray binaries, LMXBs; e.g., |Bhattacharya fc] 
van den Heuvel 1991 1 can recycle them. Owing to their 



magnetic fields suppressed by a factor ~ 10 — 10 , MSPs 
however have smaller period derivatives P oc B^/P, and 
hence longer spin-down time scales Tc ~ P/P cx (P/B)'^. 
The longer lifetimes of MSPs can therefore compensate 
for their lower birthrate so that the number of active 
Galactic MSPs can be comparable or even higher than that 



of ordinary pulsars. Lyne et al. ( 1998 1 in fact estimate 



that the number of Galactic MSPs beamed toward us 
in the radio is about equal to the number of ordinary 
pulsars, also beamed toward us, for luminosities above 
1 mjy kpc^ at 400 MHz. While the luminosity function 
of ordinary pulsars appears to flatten at luminosities 
less than about 1 mJy kpc^, these authors found no 
evidence for a flattening of the MSP luminosity function 
down to the lowest luminosities sampled of ~ 0.3 mJy 
kpc^, suggesting observationally that the underlying popu- 
lation of MSPs could in fact exceed that of ordinary pulsars. 

As we are principally interested in the high-latitude 
7— ray background from pulsars that may be confused 
with the extragalactic component, another key aspect to 
consider is the Galactic latitude distribution of the 7— ray 
emission from ordinary pulsars versus MSPs. Whereas 
young, energetic ordinary pulsars are born in and are 
concentrated close to the Galactic plane, the ages of MSPs 
generally exceed the oscillation time across the plane of the 
Galaxy (~ 100 Myr) by a large factor. Therefore, the height 
of a MSP above the Galactic plane should be uncorrelated 
with its rotational parameters. As a consequence, while the 
7— ray bright ordinary pulsars are concentrated near the 
plane, the 7— raj/ bright MSPs should have a scale height 
that is more representative of the parent population as a 
whole, and thus be more prevalent at high latitudes. This 
observation explains why the 7— ray flux distribution of 
MSPs is predicted to be close to Euclidean at the high 
end, in contrast to ordinary 7— ray pulsars which have a 
more planar distribution (Fig. [2|. The simple fact that the 
shape of the flux distribution is determined predominantly 
by geometrical considerations (see Appendix [B| allows 
us to place largely model-independent constraints on the 
contribution of MSPs to the 7— ray background (g3|. 



2.9 Model Results 

Using the Galactic pulsar population models and the 7— ray 
luminosity prescription described in §2.1| - f277| we are now 
ready to calculate the 7— ray background they produce. For 
each model, we calculate the average cumulative 7— ray pho- 
ton intensity (>100 MeV) from pulsars with absolute Galac- 
tic latitude |&| greater than a minimum value 6min: 



Ipsi{\b\ ^ fomin) = 



1 

^{\b\ ^ 6min) 



E 



(12) 



where n(|6| ^ 6min) is the solid angle subtended by latitudes 
\b\ ^ femin and the sum is over all the pulsars within that 
solid angle. Since we are interested in the high-latitude 
background that could be confused with the extragalactic 
component, we set fcniin ~ 40°. Pulsar fluxes in equation 



( 12 1 assume a 7— ray beaming fraction of unity, so that 
they spread over a solid angle 47r, consistent with the 
7— ray luminosity calibration of §2.4| Note, however, that 
this calibration was based on observed fluxes and that 
the beaming fraction factors cancel out exactly in the flux 
calculation of each pulsar. As discussed below, the 7— ray 
background however scales with the fraction of pulsars that 
have 7— ray beams encompassing our line of sight to them. 



In what follows, we will simply write Jpsr for Ip 



> 40° 



The results of our population syntheses are given in 
Table[2] For the ordinary pulsars, the model parameters are 
fixed to the ones constrained using radio observations by 
FGK06. For this model, ordinary pulsars are predicted to 
produce a small fraction, ~ 10"'^, of the high-latitude 7— ray 
background above 100 MeV. A range of parameters are 
explored for the MSP population. In ^ we will compare 
the different models with existing 7— ray observations, 
including the 7— ray background and early Fermi sources. 
In the remainder of this section, we first explore how 7psr 
depends on model parameters for the MSPs. 

Radial distribution. In our models, the Sun is at a 
distance Rq = 8.5 kpc from the Galactic center (Feast 



Ifc Whitelock| [T997| [Paczynski fc Stanek| [19981 |Reid et al 
2009). We have explored values ar = 3, 5, and 7 kpc 
(e.g., models MSPl.base, MSPl_sigr3, MSPl_sigr7) and 

for these Ipsr increases with a,., presumably because the 
bulk of the MSPs is located closer to us. 



Scale height. For the scale height of MSPs, we have 
tried values of {\z\) — 0.5, 1, 1.5, and 2 kpc (e.g., models 
MSPl.base, MSPl_z05, MSPl_zl5, MSPl_z20) and find that 
7psr increases with {\z\). This is expected, as more objects 
lie at high latitudes for larger {\z\). 

Period distribution. The 7— ray background is par- 
ticularly sensitive to the period distribution of MSPs, which 
is not surprising given the adopted scaling of the 7— ray 
luminosity of each pulsar with its period in equation Q. 
For a logarithmic slope mp,i — mp,2 = —1 of the under- 
lying period distribution, as found by |Cordes fc Chernoff| 
(19971, /par increases strongly with decreasing minimum 
period Piow as the abundance of rapidly-rotating, energetic 
MSPs increases (e.g., models MSPl.base, MSPl_Plow065 , 
MSPl_Plowl25, MSPl_Plowl5). We have also explored cases 
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Figure 2. Detailed results from our 7— ray calculations. Le/t- Map of the 7— ray sky at energies >100 MeV, its Galactic latitude profile 
at |6| ^ 40°, and the log Af — logS flux distribution for our ordinary pulsar model (FGK06). Center: Same, but for a particular MSP 
model (MSP2_base) selected for its consistency with existing observations (ijsji. Right: Same, but for a more optimistic MSP model 
(MSP2_Lrlow001). The sky maps were smoothed with a Gaussian of FWHM of 4°, similar to the EGRET point spread function. The 
solid curves in the middle panels show the best-fit disc-like (-^ 1/ sin|6|) latitude profiles. The dashed lines in the lower panels show 
different slopes /3 for the high ends of the logAf — log 5 fiux distributions, defined such that dN/dS oc S~l^. The high-latitude 7— ray 
background intensity and EGRET and Fermi source counts are summarized for all models in Table [2] 



in which the period distribution flattens to mp,i = 
at periods below a break period Pbrcak. For these cases, 
/par increases with decreasing Pbrcak at fixed Piow (e.g., 
models MSPl_Pbreakl5, MSPl J'breakS , MSPl_Pbreak5 , 
MSPl_PbreaklO), again because the abundance of short- 
period pulsars in enhanced. 

Magnetic field distribution. The magnetic field distri- 
bution is similarly important for the 7— ray background 
as the 7— ray luminosity is also assumed to be a function 
of the period derivative. In the magnetic dipole spin-down 
model, the period and period deriva tive are directly related 



models calibrated based on radio detections (^2.2). In 



to the magnetic field by P cx V PP ( Ostriker & Gunn 



1969 



[Goldreich fc Julian||1969l |Spitkovsky||2006| ). A 7-ray lumi- 
nosity law L-^ oc p'^/ ^ p--^/'' can therefore be expressed as 
Lry oc BP~^. As expected, Ip^r then increases with (logP) 
(e.g., models MSPl.base, MSPl_logB8, MSPl_logB825 , 
MSPl_logB875). At fixed (logP), it also increases with 
increasing criogs as the tail of high B pulsars becomes 
more extended (e.g., models MSPl_base, MSPl_sigB03, 
MSPl_sigB04, MSPl_sigB05). 

Radio luminosity distribution. The luminosity function 
of MSPs in the radio is also important for the cumulative 
7— ray emission predicted by the population synthesis 



models for which the radio luminosity function of MSPs 
extends to very low luminosities, the Galaxy can contain 
a large number of MSPs that are undetected in the radio, 
yet contribute to the 7— ray background. For a MSP radio 
luminosity function with logarithmic slope mL,i = — 1 
(e.g., Lyne et al.||1998 l, Jpsr thus increases with decreasing 
Pr.iow (e.g., models MSPl_base, 
Since the 7— ray back- 



low-luminosity cutoff 
MSPl_Lrlow001, MSPlXrlowlO) 
ground could in principle be arbitrarily boosted by MSPs 
which have very low radio luminosities, the background can 
be used to limit the abundance of the radio-faint sources. 
This is however not an issue for our more physical ordinary 
pulsar model, which was constrained in details in FGK06. 
Moreover, existing 7— ray observations, combined with the 
robust prediction for the shape of the MSP 7— ray flux 
distribution, constrain their contribution to the 7— ray 
background relatively independently of the model details 
(see ^ and Strong 2007 for a related approach concerning 
the Galactic diffuse emission). 

Beaming fractions. Our calculations assume a beaming 



fraction /, 



■adio,MSP 



0.5 for MSPs in the radio and f^ = l 



in the 7— rays. While the 7— ray background scales with 
the fraction of pulsars emitting 7— rays in our direction, it 
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scales inversely to the fraction whose radio beams intersect 
our sightline as a result of the calibration to the number of 
detections in radio surveys: 



to P > 16 ms and P/P < 1.25 x 10~" s"^ (|Saz Parkinson 



r _ rna 
Jpsr — Jpsr 



0.5 



(13) 



where Ip^r is the intensity obtained using the fiducial 
beaming fractions and reported in Table [2] 

Figure [2] shows detailed results from our 7— ray calcu- 
lations. The left column shows a map of the 7— ray sky, its 
Galactic latitude profile at |6j ^ 40°, and the log TV — log 5 
flux distribution for our ordinary pulsar model. The solid 
curve in the middle panel shows the best-fit disc-like 
(~ l/sin|6|) latitude profile. The central column shows 
the same quantities, but for a particular MSP model 
(MSP2_base) selected for its consistency with existing 
observations, to be discussed in the next section. The right 
column also shows the same quantities but for a more 
optimistic MSP model (MSP2_Lrlow00l). For the sky maps, 
sources with fiux S > 10~^ ph s^^ cm~^ (the approximate 



point-source sensitivity of EGRET; see were omitted 



to avoid bright spots corresponding to resolved sources. 
The high-latitude 7— ray background intensity is given in 
Table [2] for all models explored. 



3 COMPARISON WITH EXISTING 
GAMMA-RAY OBSERVATIONS 

Next, we explore the constraints that existing 7— ray ob- 
servations put on the contribution of pulsars to the 7— ray 
background. 



3.1 Existing Detections 

An important test of our 7— ray pulsar population models is 
provided by the number of detections by EGRET and early 
Fermi observations. Specifically, the number of high 7— ray 
flux pulsars predicted by a given model must be consistent 
with existing EGRET and Fermi results for the model to 
be viable. The strongest constraints are provided by the 
first three months Fermi catalog (lAbdo et al. 2009c I, in 



which a larger fraction of the sources are identified than in 
the EGRET data, and so we focus on these. 

Unfortunately, the comparison is complicated by the 
complex selection effects in detecting 7— ray pulsars. For 
example, three different approaches to detecting pulsars 
are adopted by Fermi (Saz Parkinson 20091. First, 7— ray 



pulsations are searched for toward pulsars that are already 
known, usually in the radio. Because 7— ray sources are 
inherently dim, it is easier to detect pulsations when their 
period is known from measurements at other wavelengths. 
Second, blind searches for 7— ray pulsars are carried out 
using a time-differencing technique toward unidentified 
EGRET and Fermi sources. These searches are important 
since many pulsars appear to be radio-quiet (e.g., the 
Geminga and CTA 1 pulsars) and may only be detectable 
in the 7— rays. They are however generally less sensitive 
when the period is not known a priori and currently limited 



f2009l). Third, many of the blind searches actually target 
well-studied astrophysical sources (including EGRET 
unidentified sources, supernova remnants, and TeV sources 
discovered by ground-based Cerenkov telescopes) which, 
to various degrees of confidence, are suspected to contain 
pulsars. The detection threshold is also a function of 
position on the sky, with the diffuse background being 
stronger near the Galactic plane (e.g., Abdo et al.|2009c l. 



To avoid the large uncertainties from the complicated 
selection function, we therefore base our discussion on 
the flux complete subset of the analysis of the \b\ 10° 



sources reported by Abdo et al. (2009b I. This corresponds 
to a flux threshold of S > 1.25 x 10"'^ ph s"^ cm"^ at 
energies >100 MeV. We denote the number of simulated 
sources with flux above this threshold and with latitude 



|b| ^ 10° by iV„ 



^ 10°). Among the sources with flux 



above this threshold, two are new pulsars discovered by 
Fermt (LAT PSR J1836+5925 and LAT PSR J0007-h7303) 
and two more are unidentified. While 40 are associated at 
high-confidence with AGNs, the high-confidence associa- 
tions are only defined as having probability >90% of being 
genuine, so that a few could be spurious. In addition, PSR 
J0357. 5-1-3205 is a previously-known MSP with 7— ray flux 
just below the completeness threshold, but consistent with 
a value above the threshold within the uncertainty. Because 
the 7— ray beams of pulsars appear generally wider than 
their radio beams and because many pulsars have not yet 
been detected at other wavelengths, we can expect at least 
as many 7— ray pulsars to be unidentified as sources with 
known counterparts. It is thus likely that the sources in the 
flux limited sample that are not associated with a known 
AGN are pulsars. Conservatively and accounting for the 
fact that some of the potential sources may not be MSPs, 



we fiducially take the allowed simulated 



> 10° 



to be 4, with a significant statistical uncertainty. At the 2a 
level, 8 sources are allowed. The predicted A^compd&l ^ 10°) 
for each model is given in Table [2] 



3.2 Background Intensity and Latitude Profile 

The extragalactic 7— ray background is usually measured by 
subtracting from the total intensity a Galactic foreground 
model accounting for 7— rays produced by cosmic-ray 
nucleons interacting with nucleons in the interstellar 
gas, bremsstrahlung by cosmic-ray electrons, and inverse 
Compton interaction of cosmic-ray electrons with ambient 
low-energy interstellar photons (e.g., Sreekumar et al.|1998 



[Strong et al.|[2004| . A potential pulsar contribution would 
therefore be lumped into the extragalactic component Ix- 
Existing extragalactic background estimates can therefore 
be used to constrain the Galactic pulsar population and its 
7— ray emission from the requirement that their integrated 
7— ray intensity must not exceed these. Table[l]lists existing 
estimates of the extragalactic 7— ray background above 100 
MeV, along with the corresponding spectral shape in the 
cases where it was also measured. 

The original analysis of the final EGRET data by 
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Table 1. Extragalactic 7-ray (>100 MeV) Background Estimates 



Reference 


Ix 


Spectrum 


Instrument 




10~^ ph s~^ cm~^ sr~^ 







Th ompson fc Fich tel 1 1982 I 
JOsborne et al' (wd^ 
Sreekumar et al. ( 199 
;Keshet et al (2004]l 
jStrong et aL| | |2004| 



1.3 ±0.5 

1.10 ±0.05 
1.45 ±0.05 

< 0.5 (99% conf.)b 

1.11 ± 0.01 



: 2.35 



-0.4 



P - — _0.3 

p = 2.11 ± 0.05=^ 
p = 2.10 ± 0.03^ 

positive curvature'^ 



SAS-2 
EGRET 
EGRET 
EGRET 
EGRET 



"Power-law photon index. 

''These authors found that the 7-ray background Galactic latitude profile is disc-like (~ 1/ sin|6|) for 40° ^ |b| ^ 70° and even steeper 



for |6| > 70°, with an average intensity at the Galactic poles > 86°) /pole = 

foreground subtraction 



1.20 ± 0.08 X 10-5 ph s-i cm"^ sr'^ 
Not well described by a power law. The positive curvature may be affected by systematic uncertainties. 



before Galactic 



ph s 



cm ^ sr 



Sreekumar et al. ( 1998 1 found Ix = (1.45 ± 0.05) x 10" 



while a more recent analysis by Strong 



et al. (2004) using an updated foreground model found a 
lower value Ix = (1.11 ±0.01) x ID"'' ph s"^ cm'^ sr'^ At 
present about 7% of the extragalactic background has been 



resolved into blazars by early Fermi observations (Abdo 
et al. 2009b'), though this fraction is likely to increase as 



deeper integrations are obtained. An upper limit to the 
pulsar contribution is thus /psr 1 x 10^^ ph s~^ cm^^ sr^^. 

In a more model-independent analysis based on corre- 
lations with Galactic tracers and on the Galactic latitude 



profile of the 7— ray background, Keshet et al. ( 2004 1 find 
Ix < 0.5 x 10^5 ph s"^ cm~^ sr^ 
level, with evidence for an even 



at the 99% confidence 
lower intensity. The 



the measurements of Sreekumar et al 



et al. 



discrepancy between the Keshet et al. ( 2004 1 analysis and 

(11998!) and [Strong 



( 2004 1 could be understood if pulsars were picked 
up as a Galactic component in the [Keshet et al.| ( [2004[ ) 
analysis. The 7— ray Galactic latitude profiles predicted for 
our pulsars models (Figure [2|, with their nearly disc-like 
behavior ~ 1/ sin|6|, are in fact tantalizingly similar to the 
one measured at |fo| 40° (see Figure 3 of Keshet et al. 



2004) 



3.3 Background Spectrum 

The sources that dominate the isotropic component of the 
7— ray background must have spectra that are consistent 
with the latter. While there remains significant variance 
between the different estimates (Table [T]), these indicate 
a spectrum with approximately equal power per energy 
decade (a power-law photon index p « 2) , or slightly steeper 
and possibly with positive curvature, between 100 MeV and 
100 GeV. In fact, the similar spectra of blazars were used 
as an argument in favor of attributing most of the isotropic 
component to these sources ( [Sreekumar et aL][1998[ [Strong[ 
et al.|2004| ). 



Are the 7— ray spectra of pulsars also compatible with them 
making a significant contribution to the high-latitude 7— ray 
background? The relatively fiat spectra of pulsars, up to a 
certain energy, are indeed in agreement with the EGRET 
7— ray background spectrum measurements. However, of the 



few pulsars with 7— ray spectra, many exhibit high-energy 
cutoffs around a few GeV. For instance, the spectrum of 
the Vela pulsar is exponentially suppressed at around 2.9 



GeV (Abdo et al. 2009a). The first MSP with a spectrum 



measured by Fermi, PSR J0030+045, also exhibits a cut off 
at around 1.7 GeV ( Abdo et al.|2009b ). Such a pronounced 
cutoff is not observed in the 7— ray background spectrum 
measured by [Sreekumar et al. |l998j ) at these energies 
and would place a strong constraint on the contribution 
of pulsars if it were typical of the pulsars dominating the 
background component. Fermi measurements of the 7— ray 
spectra of more MSPs will therefore be extremely valuable 
in quantifying their contribution to the 7— ray background. 

Note, however, that in their analysis of the same EGRET 
data with an improved Galactic cosmic-ray foreground 
model Strong et al. ( |2004 1 measure a residual 7— ray 
background spectrum that is inconsistent with both a 
simple power law and with the original [Sreekumar et aL] 
(|1998|) analysis. The Strong et al. ( 2004 ) measurement is in 



fact indicative of an upturn in the spectrum around a few 
GeV. It is thus possible that pulsars contribute a significant 
fraction of the background up to energies of a few GeV even 
if their spectra are suppressed at these energies, if other 
sources provide the photons observed at higher energies. It 
will be necessary to also obtain a firmer measurement of 
the background spectrum before stronger constraints can 
be put on the contribution of pulsars to this background 
from measurements of their individual spectra. 



A preliminary analysis of the Fermi data 
presentation by M. Ackermann at 



the 



[http://www- 

conf. slac ■ Stanford . edu/t evpa09 / [ ) indicates an isotropic 



7— ray background spectrum that is steeper than the 



original EGRET analysis by [Sreekumar et al. ( 1998 ) with 
p « 2.45 between 200 MeV and 50 GeV, and of lower 
intensity, in better agreement with the [Strong et al.[ ( [2004[ ) 
normalization. As the error bars converge, the presence 
or absence of a feature in this spectrum around few GeV 
will strongly constrain the contribution of MSPs if their 
spectra are generically suppressed at these energies. The 
preliminary Fermi error bars at a few GeV, where the 
observed pulsar spectra peak relative to the background 
spectrum, are large enough to comfortably accommodate 
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a 5 — 15% contribution from MSPs, as allowed by other 



constraints {[ 3.4 1 



3.4 Which Models Pass the Tests? 

Among the models explored (Table [2|, which ones are 
compatible with the existing data? The basic requirements 
are: 1) iVcomp(|&l > 10°) < 4 and 2) Ip^, < 1 x 10"^ ph s^^ 



With a fractional high-latitude 7— ray background contri- 
bution of 10""^, our model for ordinary pulsars appears 
quite consistent with the data. Because the 7— ray bright or- 
dinary pulsars are concentrated near the Galactic plane, the 
|6| > 10° data provide weak constraints. We have checked 
that the model predicts only a total of 17 sources on the 
sky above a flux threshold of 1 x 10~^ ph s~^ cm~^, roughly 
the average EGRET point source sensitivitjj^ ( [Gehrels 
& Michelson 19991. This is in good agreement with the 



six pulsars detected by EGRET and the 170 unidentified 



sources in the third catalog (Hartman et al. 19991. Most 



of our MSP models can however already be excluded on 
the basis that they overproduce either Ncoinp{\b\ ^ 10°) 
or the 7— ray background. There are however viable and 
interesting models. Model MSP2_base (with Piow = 1.5 ms 
and {\ogB) = 8.0), whose detailed results are plotted in 
Figure [2j predicts Ncomp{\b\ ^ 10°) = 4 and a high-latitude 
7— ray background from MSPs Jpsr ~ 0.08 x 10~^ ph s~^ 
cm~^ sr~^. This corresponds to a sizable fraction of the 
extragalactic value: slightly less than ~ 10% if the mea- 
surements obtained by subtracting a cosmic ray model are 
correct, but ~ 20% or more if the more model-independent 
analysis of Keshet et al. (20041 is more accurate (Table [T]). 
It is interesting that the MSPs that are the most luminous 
in the 7— rays are likely to be those with the highest E 
and therefore the most likely to ablate their companions 
and be eclipsed by the vaporized material. Such pulsars, 
which would be difficult to identify in the radio or in 
blind 7— ray searches, could remain 7— ray sources either 
by transmitting most of their intrinsic 7— ray luminosity 
or by secondary 7— rays produced in the interaction with 
the obscuring material ( |Tavani||l991[ ). The more optimistic 
MSP model MSP2_Lrlow001, whose detailed results are 
shown in the right column of Figure [2] overpredicts the 
number early Fermi detections, but illustrates how MSPs 
could potentially contribute an even large fraction of the 
background under reasonable assumptions. 

The fact that many of our MSP models, with quite 
reasonable population parameters based on the litera- 
ture, overproduce the high-latitude 7— ray background is 
particularly notable, as it implies that the background 
puts interesting constraints on the global Galactic MSP 
population. We currently only know of ~ 100 MSPs in the 
Galactic field, although there could be ~ 10^ or more in 



total (Lyne et al. 19981, many of which would be unde- 



tectable in the radio because of beaming. Although this 



^ Near the Galactic plane, the point source detection threshold 
is somewhat higher owing to the background, so this estimate is 
optimistic. 



may suggest that our MSP 7— ray background predictions 
are wildly uncertain, this is not the case. In general, 
the number of resolved point sources above a given flux 
threshold is not independent of the integrated flux from 
the parent population. In Appendix |A] we derive this 
relation and analytically estimate the contribution of MSPs 
to the high-latitude 7— ray background from the number 
of resolved point sources in the early complete Fermi 
sample. We find Jpsr ~ 0.125 x lO"'^ ph s"^ cm"^ sr"^ 
for A^comp(|fe| 10°) = 4, in good agreement with our 
numerical calculations for models with consistent numbers 
of resolved sources. This estimate relies on the slope j3 of the 
log A'^ — log S* distribution (defined such that dN/dS oc S~'^), 
which is empirically found to be nearly Euclidean (/3 = 2.5) 
at the high-fiux end (Figure [2|. In Appendix [b] we show 
that this slope should be robust to the details of the 
population models, as we have empirically verified to be 
the case for all of our MSP models. It is also worth noting 
that Story et al. (20071, who performed more extensive 



calibrations to radio data for MSPs and treated beaming 
geometries in great detail, predict a nearly identical shape 
for the 7— ray flux distribution (see their Fig. 10). Note 
again that both the small-number statistical uncertainty on 
Acompd&l ^ 10°) and the possibility of a higher value owing 
to spurious associations of some sources with AGNs might 
allow an even more significant contribution of MSPs to the 
7— ray background. 



4 IMPLICATIONS FOR PULSAR PHYSICS 

The previous section was chiefiy concerned with the con- 
straints that can be put on the pulsar population as whole, 
particularly MSPs, using the 7— ray background. These also 
have implications for the physical processes that govern pul- 
sars. 



4.1 Maximum Spin Frequency 

If 7— ray emission from pulsars, including MSPs, is in 
fact a function of the period and period derivative of the 
pulsar of the type in equation ([3|, then the measured 
7— ray background puts a lower limit on the maximum 
spin frequency (or minimum spin-period Piow), as argued 
!|3] This has interesting consequences arising from the 
physics that limit this spin frequency: the break-up spin 
rate, accretion spin-magnetic equilibrium, and gravitational 
wave emission. We discuss these below. 

A strict upper limit on the spin frequency of neutron 
stars is the angular velocity at which centrifugal accelera- 
tion would disrupt the star, usually called the "break-up 
limit." This break-up frequency depends on the equation of 
state of the ultradense matter at the cores of neutron stars. 
A firm upper limit of « 3 kHz is set by the requirement 
that the sound speed within the star cannot exceed the 
speed of light. More realistic equations of state predict 
lower maximum spin rates in the 1500—2000 Hz range 



(Cook et al. 1994 Haensel et al. 1999 Lattimer & Prakash 



2001 Chakrabartyl2008 1 . In order to accommodate a radius 



of 8—12 km for a 


1.4 Mr neutron star (e.g., iThorsett & 


Chakrabarty|l999 


6zel|2006 


Guver et al.|2008 


Ozel et al. 
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20091, a break-up limit greater than about 1500 Hz is 



needed. 

It took 25 years from the discovery of the first MSP, 
B1937-I-21 with a spin period of 1.6 ms, before a more- 
rapidly rotating pulsar, J1748— 2446ad with a period of 1.4 
ms, was found (Hessels et al. 20061. While 'Hcssels et al 



( 2006 1 argue that the difficulty of detecting such a pulsar 
suggests that more with comparable or shorter periods 
may exist in the Galaxy, observations of nuclear-powered 
millisecond X-ray pulsars, which during outbursts are 
bright enough to be detected throughout the Galaxy and 
are not affected by the same selection effects, however 
indicate that 1.4 ms is close the minimum spin period 
(e.g. 



like the ones we have explored could be used to better 
constrain the short-end of the period distribution using the 
7— background and thus help determine the physics that 
limits pulsar spin rates. 

4.2 Gamma-Ray Emission Mechanisms 

The most obvious question, of course, is how can the 7— ray 
background inform us on the pulsar 7— ray emission mech- 
anisms? Two main classes of 7— ray emission mechanisms, 
reviewed by Harding et al. (20071, have been proposed. 



Polar cap models assume that particle acceleration takes 
place near the neutron star magnetic poles and that 7— ray 
emission results from cascades initiated by magnetic pair 



Chakrabartyl |2008| . The fact that this putative production in the strong magnetic field (e.g. [Daugherty 



minimum spin period is a factor of two or more above the 
break-up limit for favored equations of state is interest- 
ing as it may be evidence that a different mechanism is 
responsible for the dearth of pulsars spinning at higher rates. 

The most conservative explanation for a maximum 
spin frequency above the break-up limit is rooted in the 
origin of MSPs in LMXBs. In this recycling scenario, the 
binary companion to the neutron star fills its Roche lobe 
and transfer of orbital angular momentum to spin angular 
momentum ensues from accretion of material onto the 
compact object (e.g 



Bhattacharya & van den Heuvel 



1991). Because the neutron star is strongly magnetized. 



the Alfven radius Ra, at which the kinetic energy density 
of the accreting material becomes comparable to the local 
magnetic energy density, is well outside of the neutron star 
surface. As a result, the accretion flow exerts viscous torques 
at Ra until an equilibrium spin period, approximated by 
the angular Keplerian velocity at Ra, is reached: 



Peq ~ 1.9 ms 



B 



109 G 



6/7 



M 

MBd 



-3/7 



R* 



10 km 



16/7 



(14) 

where M / M^d is the mass accretion rate as a fraction of 
the Eddington limit and Rt, is the radius of the neutron 
star (e.g., Ghosh fc Lamb|1992 l. A short-period cutoff that 
is the result of magnetic spin equilibrium would therefore 
constrain the joint distribution of B, M, and R. A sharp 
period cutoff seems unlikely in this scenario as, for example, 
there is no obvious reason for B to be strictly limited to a 
particular range; rather, the P — P diagram is suggestive of 
a broad B distribution. 

Another interesting possibility is that accretion spin- 
up torques may be mitigated by angular momentum losses 
owing to gravitational radiation. These could occur via 
a number of proposed mechanisms, including r— modes 
or a quadrupole moment induced either by accretion or 
a toroidal magnetic fleld (e.g.. Wagoner] 



1984 



Andersson 



eraL][T999l [Ushomirsky et al.||2000| |Cutler||2002| [Melatos 



& Payne 20051. An attractive aspect of gravitational wave 



braking is that the spin-down torque scales with the fifth 
power of the spin frequency, oc lo^ , and would thus more 
naturally produce a sharp suppression at short spin periods. 

If Fermi observations can provide a reliable 7— ray 
emission law for MSPs, then population synthesis models 



& Harding 1982 19961. Outer gap models instead assume 



that the particle acceleration takes place in vacuum gaps 
that form in the outer magnetosphere (e.g., Cheng et al. 
1986 Romani|1996 1. The slot gap model is a variant of the 
polar cap model in which acceleration occurs in a 'slot gap', 
a narrow region bordering the last open field line in which 
the electric field is unscreened. The high-energy emission 
beam of the slot gap cascade is a hollow cone with a much 
larger opening angle than that of the polar cap cascade 
emission. 

7— ray observations provide a few handles on the pul- 
sar emission mechanisms. While our empirical approach 
to modelling 7— ray emission does not allow us to make 
quantitative statements, we briefiy outline discriminating 
statistics that further studies within the framework of 
particular models will be able to make precise (for existing 
work in this direction, see|Gonthier et al.|2002l |2004[ |Story| 



|et al.|2007| [Harding et al.|2007 l 



7— ray beaming fraction. Outer gap models predict 
wide 7— ray beams, with beaming fractions of order unity, 
while polar cap models predict much smaller beams com- 
parable to those in the radio. Therefore, outer gap models 
predict a large ratio of radio-quiet to radio-loud 7— ray 
pulsars and a stronger integrated 7— ray background. Slot 
gap models are intermediate. 

7— ray spectrum. The different classes of models predict 
different high-energy cutoffs. Polar cap models, for example, 
can have hyper-exponential cutoffs that are a signature of 
7 — 5 pair attenuation of low-altitude emission. If MSPs 
are found to contribute a significant fraction of the 7— ray 
background, then the spectrum of the latter (possibly after 
subtracting the blazar component) will inform us about 
the average MSP spectrum. This information would be 
complementary to observations of individual MSPs, since 
it is conceivable that the brightest MSPs have different 
properties from the bulk of the population. 

Dependence on magnetic and rotational parame- 
ters. Whether MSPs do contribute a significant fraction 
of the 7— ray background depends in part on whether the 
empirically-motivated relation in equation (|3| holds for 
all pulsars. Given a MSP population model and if a large 
fraction of the high-latitude 7— ray background is found to 
arise from blazars (or other non-pulsar sources), one can 



14 Claude-Andre Faucher-Giguere and Abraham Loeb 



Ordinary Pulsars (FGK06) Viable Model IVISP2_base Optimistic IVIodel MSP2_Lrlow001 

>100 MeV, S]„„ = (10 degf >W0 MeV, n^,, = {W degf >100 MeV, n^, = (W degf 




Figure 3. Illustration of how the pixel count PDF can distinguish between different sources of the 7— ray background. We show the 
normalized PDF for pixels at latitudes |b| ^ 40° for an exposure X = 2.5 X 10^^ cm^ s, corresponding to a 5-year all-sky Fermi survey. 
We consider two pixel size and energy range combinations: >100 MeV, fipix = (10°)^ and >1 GeV, fipix = (1°)^- The left column shows 
the results for our ordinary pulsar model, the central column the results for a particular MSP model (model MSP2_base) selected for its 
consistency with existing observations (Sjsj), and the right column the results for a more optimistic MSP model (model MSP2_Lrlow001). In 
each panel, the dashed curve shows the Poisson distribution expected if the background is truly diffuse and uniform on the sky, with the 
same total number of counts. Pulsar models are clearly distinct, with much broader PDFs. The pixel count PDF could also be computed 
for any particular blazar model and compared with pulsar predictions. 



constrain the range of applicability of any particular 7— ray 
luminosity law. 



5 FLUCTUATIONS IN THE 7 RAY SKY 

Most of this work has focused on average 7— ray back- 
ground observables, such as its intensity and spectrum. 
Fluctuations in the 7— ray sky can however also be fruitfully 
exploited to discriminate between different origins of the 
7— ray background, even if its sources are unresolved. 
In this section, we show how a particular statistic, the 
probabihty distribution function (PDF) of the number 
of 7— ray counts in sky pixels, can be used to distin- 
guish the 7— ray background of pulsars from truly diffuse 
emission and potentially also from other sources like blazars. 

We imagine dividing a portion of the sky of interest 
into A^pix angular pixels each of solid angle fipix. The region 
of the sky is then observed for an exposure X in units of cm'^ 
s, reflecting a combination of effective area and integration 
time. At the end of the observation period, each pixel has 
collected a number of 7— ray photons A'^^ . By combining the 
data from all the pixels, we can estimate P{Nj), the PDF 
of the number of photons collected by an individual pixel. 
The method takes advantage of the ability of the telescope 
to measure the direction of arrival of the 7— ray photons. 
Fermi, for instance, has an energy-dependent point spread 
function, with an angular resolution that improves with 



increasing energy. We may thus potentially gain with this 
approach by considering only photons above a given energy 
threshold. The 68% containment angle for Fermi varies 
from approximately 5° to 0.5° from 100 MeV to 1 GeVQ 

5.1 Analytical Derivation 

We first assume that the 7— ray background originates from 
a population of uncorrelated discrete sources, with a photon 
flux distribution dN/dS above a prescribed energy thresh- 
old, and a total of Nsrc sources in the sky area considered. 
By the law of total probability, 
00 

Pi^i) = ^(^■y\ ^ pix)P(fc src in pix), (15) 

where 'fc src in pix' represents the case that exactly k sources 
lie in the pixel. The distribution of the number of sources 
within each pixel is a simple binomial: 

P{k src in pix) = j Ppix(l - Ppix)^'""'', (16) 

where ppix = is the probability that a source lies within 
any particular pixel. If a pixel contains a single source of 
photon flux 5, the number of photons collected after an 

http://www-glast.slac.stanford.edu/software/IS/ 
glast _lat _performance . htm 
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exposure X has a Poisson distribution with a mean equal 
to the expected number at the end of the exposure, SX. 
Taking into account the photon flux distribution, 



P{N.y \ 1 src in pix) — 



iVsrc 



^^dNe-^-jSXr^ ^ (17) 



dS 



7V-, 



Note that because of the integration over the photon flux 
distribution, this PDF in general differs strongly from a 
Poisson distribution. 

If fc > 1 sources lie in the pixel, A'^-, is distributed 
like the sum of independent random variables with the 



distribution in equation (171. -P(A^7| k src in pix) can 



then be expressed as a convolution of PDFs of the form 
of the one-source case. For small k, the PDF will also be 
non-Poissonian. For large k, however, the central limit 
theorem implies that the distribution will tend to become 
Gaussian and therefore lose discriminating power. A good 
choice for the pixel size is therefore one for which (given 
the instrumental limitations) each pixel contains a small 
number of bright sources, so that P{N.y) is determined by 
the small-fc P{N-,\ k src in pix) terms in equation (15 I. 



For comparison, consider the case of a 7— ray back- 
ground with exactly the same average intensity and 
spectrum as the discrete case above, but completely diffuse 
and uniform. Then P{Nj) is a simple Poisson distribution 
with mean equal to the expected number of photons per 
pixel at the end of the exposure, (N^) — {Nsic/Npi,^)(S)X . 
Therefore, even for same intensity and spectrum, the pixel 
count PDF can distinguish between between a discrete 
and diffuse origin for the 7— ray background if the angular 
resolution of the instrument is sufficient. This idea can be 
extended to distinguish competing models that consist of 



discrete sources. As equation (171 shows, P{N^) depends 
on the full photon flux distribution (and total number) 
of the sources. The PDF could therefore help distinguish, 
for example, pulsar and blazar models, or even between 
different pulsar models. 

5.2 Numerical Models and Application to Fermi 

We may use the numerical models of S|2]to compute the pixel 
count PDF expected from pulsars. We focus on the |6| ^ 40° 
pixels of a 5-year Fermi all-sky survey, which corresponds 
to an exposure of about X = 2.5 x lO^'^ cm'^ s (Baltz et al. 



2008). We explore two combinations of pixel size and energy 
range resolvable by Fermi: >100 MeV, fipix = (10°)^ and 
>1 GeV, fipix = (1°)^- For these numerical calculations, we 
divide the sky in our simulations into pixels and simulate 
an actual observation. The results are shown in Figure |3j 
with our ordinary pulsar model on the left, one of our viable 
MSP models (model MSP2_base) at the center, and a more 
optimistic MSP model (model MSP2_Lrlow00l) on the right. 
In each panel, the dashed curve shows the Poisson distribu- 
tion expected if the background is truly diffuse and uniform 
on the sky (as may be expected, for example, in some dark 
matter annihilation scenarios), with the same total num- 
ber of counts. Pulsar models are clearly distinct, with much 
broader PDFs. The pixel count PDF could also be computed 
for any particular blazar model and compared with pulsar 
predictions to help identify which are the dominant sources. 



Note, however, that different population models may yield 
different PDFs and so a proper analysis of the data will need 
to explore a range of both pulsar and blazar models in order 
to robustly distinguish between the two hypotheses. 



6 CONCLUSIONS 

We have estimated the contribution of Galactic pulsars to 
the diffuse 7— ray background, focusing on the high-latitude 
part (|6| ^ 40°) that could be confused with an extragalac- 
tic component. We have considered both the population of 
ordinary pulsars and the population of recycled MSPs. For 
the ordinary pulsars, the underlying population parameters 
were constrained in detail by the radio population study of 
FGK06 and taken as fixed. Since the MSP population is 
not as well constrained observationally due to more severe 
and complex selection effects, we instead explored a range 
of models with fiducial parameters based on the literature, 
also calibrated to the number of detections in radio surveys. 
We examined how the predicted 7— ray background, as well 
as the predicted number early Fermi detections, depend 
on the population parameters as a guide to how 7— ray 
measurements can usefully constrain the MSP population. 
To assign 7— ray luminosities to synthetic pulsars, we 
adopted a simple empirical prescription motivated by 
EGRET results and which also provides a good fit to 
early Fermi and AGILE observations, including a some 
MSPs. According to this prescription (eq. |3| that upcoming 
Fermi observations will help refine, the 7— ray luminosity 
of each pulsar is a simple function of its period and period 
derivative, up to a maximum fraction of the total power 
available from spin-down. 

We have found that while ordinary pulsars are ex- 



pected to contribute a relatively small fraction 



10" 



of the high-latitude 7— ray background, MSPs could in 
principle easily contribute a significant fraction of the 
measured high-latitude intensity, Jx ~ 1 x 10"^ ph s~^ 
cm~'^ sr~^. It is therefore possible to use measurements of 
the 7— ray background to usefully constrain this population 
(^; many population models are in fact already ruled out 
by existing 7— ray observations. Models that are consistent 
with early Fermi results suggest that MSPs could con- 
tribute ~ 5 — 15% of the high-latitude 7— ray background 
if the measurements obtained by subtracting a cosmic ray 
model are correct, but up to ~ 20% or more if the more 



model-independent analysis of Keshet et al. ( 2004 1 is more 



accurate. This conclusion is actually independent of the 
assumed details of the MSP population and rests principally 
on the simple fact that the ages of MSPs generally exceed 
the oscillation time across the plane of the Galaxy (~ 100 
Myr) by a large factor. Unlike the 7— ray bright ordinary 
pulsars which are concentrated near the plane, the 7— ray 
bright MSPs are equally prevalent at all latitudes (modulo 
the population scale height). This results in a logA^ — \ogS 
fiux distribution that is nearly Euclidean at the high-end 
for MSPs, which can be used to relate the integrated 
background to the number of resolved point sources above 
a given fiux threshold. 

For our adopted 7— ray prescription, the 7— ray back- 
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ground intensity is particularly sensitive to the period of 
the most rapidly-rotating MSPs. If this maximum spin rate 
can be constrained to be significantly below the break-up 
upper limit, it would provide information on the processes 
that limit the spin frequency. This would constrain the 
conditions of spin-magnetic equilibrium in which MSPs 
are thought to acquire their high rotation frequencies, 
and may be indicative of braking by gravitational wave 
emission, which would have important ramifications for 
direct detection efforts by Z//GO[^ 

Given that pulsars, and MSPs in particular, appear 
quite capable of contributing significantly to the 7— ray 
background, it will clearly be worthwhile for upcoming 
Fermi studies to check and refine the 7— ray prescription 
we have used in our calculations. More detailed population 
synthesis studies of MSPs in the radio, for which there 
has been relatively little work, will also be very useful. 
If the 7— ray contribution of Galactic MSPs proves to be 
important, it will be necessary to also subtract it in order to 
obtain an accurate measurement of the extragalactic back- 
ground. MSPs will have to be considered when attempting 
to detect a putative component from dark matter annihila- 
tion, which is likely to produce a subdominant signal. As 
we showed in ^ fluctuations in the 7— sky (which depend 
on the full log TV — log S flux distribution) can be fruitfully 
exploited to distinguish between different sources of the 
7— ray background. As the telescope exposure increases, 
more of the unresolved discrete sources contributing to the 
fluctuations will be individually detected. 

As a flnal recommendation for future observations, we 
emphasize the importance of blind searches since a substan- 
tial fraction of the 7— ray background could originate from 
pulsars that are not visible in the radio. This is especially 
true for MSPs that can be difficult to detect in the radio 
(either because of instrumental limitations, dispersion or 
scattering by intervening interstellar electrons, or because 
of orbital acceleration or eclipses in binary systems) but 
may be more prevalent as 7— ray sources. At present, blind 
time-differencing searches for periodicity by Fermi are 
limited to P > 16 ms and P/P < 1.25 x 10"" s'^ |Sm 



Parkinson||2009 ). If possible (given the limitations imposed 
by the small numbers of counts) the Fermi search space 
should be increased to allow the detection of shorter-period 
MSPs and also for acceleration in binary systems, as 
these should be both numerous and may be significant 
contributors to the diffuse background. 
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APPENDIX A: RELATION BETWEEN 
RESOLVED SOURCES AND THE 
INTEGRATED BACKGROUND 

The contribution of a source population to the integrated 
background intensity, I, is in general not independent of the 
number of resolved sources above a given flux threshold, 
N{> Sth)- If the log A'' — logS flux distribution is known, 
the two are in fact a direct function of one another. While 
the actual full flux distribution of the Galactic pulsar 
population is not yet empirically measured, we can combine 
the theoretical predictions of !|2] with early Fermi results 
to more model-independently quantify the contribution of 
pulsars to the 7— ray background intensity. 

For simplicity, let us assume that the pulsar 7— ray 
flux distribution has a simple power-law form above a lower 

limit S'min, 



dN 

Is 



AS- 



if iS* ^ 5'miii 

otherwise 



Then, provided that Sth > 5'miii and /? > 2, 

roc A 

Ni> Sth) = / dN= j^Sl-", 
while the integrated flux from all the sources 



Stot — I 

Jo 



dNS - 



(3-2 



2-/3 



(Al) 



(A2) 



(A3) 



Deflning the average background intensity over a solid angle 
Q as Ipsr = Stot/^, we obtain the simple relation 



1-/3 



N{> 5th). (A4) 



To compare with existing observations, we can consider the 
portion of the sky at Galactic latitudes |6| ^ fcmin, for which 
Q — 4-k{1 — sinfeniin)- We base our discussion on the analysis 
of the l&l 10° sources detected by Fermi in its first three 



months of operation reported by Abdo et al. ( 2009bJ (see 
also 1 3.4 1. This sample is complete for fluxes S > 1.25 x 10~^ 
ph s~ cm~^ at energies >100 MeV. Among the sources 
with flux above this threshold, two are new pulsars dis- 
covered by Fermi (LAT PSR J1836-)-5925 and LAT PSR 
J00074-7303) and two more are unidentifled. While 40 are as- 
sociated at high-confidence with AGNs, the high-confidence 
associations are only defined as having probability >90% 
of being genuine, so that a few could be spurious. In addi- 
tion, PSR J0357.5-f3205 is a previously-known MSP with 
7— ray flux just below the completeness threshold, but con- 
sistent with a value above the threshold within the uncer- 
tainty. Because the 7— ray beams of pulsars appear gener- 
ally wider than their radio beams and because many pulsars 
have not yet been detected at other wavelengths, we can ex- 
pect at least as many 7— ray pulsars to be unidentifled as 
sources with known counterparts. It is thus likely that the 
sources in the flux limited sample that are not associated 
with a known AGN are pulsars. Conservatively and account- 
ing for the fact that some of the potential sources may not 
be MSPs, we flducially assume that there are four 7— ray 
MSPs in the flux-limited sample. For an Euclidean source 
distribution (/3 = 2.5; see Appendix [B] below) , the Fermi 
resolved sources thus allow for a contribution of pulsars to 
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the |6| ^ 10° 7— ray background at energies > 100 MeV of 



1 -2 -1 



> 10 ) ~ 0.5 X 10"" ph s"" cm"^ sr 



Sn 



10-10 ph s-i cm-2 



Sth 



1.25 X 10-'^ ph s-i cm-2 
fN{> Sth) 



V 4 



(A5) 



with a 50% uncertainty from small-number statistics only. 
In the above estimate, the 5min factor was approximately 
set to the turn over of the MSP log A'^ — log S distribution 
seen in our MSP population synthesis models (*|2|. If the 
Fermi flux-limited sample contains spurious associations of 
sources with AGNs because of positional coincidence, then 
the true 7V(> Sth) for MSPs could be higher and we would 
therefore expect a larger contribution from these sources to 
the background intensity. 



In this work, as in the study of [Keshet et al.| ( |2004[ ), 
we defined the intensity of the high-latitude 7— ray 
background as the |&| ^ 40° average {[2.9\. However, 



the Fermi fiux limited sample on which we base our 
point source comparison covers |fe| 10° (Abdo et al. 
2009b I. Assuming that the latitude profile is disc-like 



(~ l/sin|6|; Figure [2|, there is a simple expression for the 
ratio of the average intensities at |fo| ^ &min for different 
values of femin- This expression follows from the result 
{(sin > 6min) = In [(sin |6min|)^''"] for the solid 

angle average: 



(A6) 



I{\b\ > &min,i) ^ In [(sin|bn,in,i|)-i] 

Ii\b\ > &min,2) ln[(sin|&n,in,2|)-l] 

Numerically, I{\b\ > 40°) ^ 0.25/(|6| > 10°) and the numer- 
ical factor in the analytic estimate of equation ( A5 1 becomes 
0.125 X 10-=^ ph s"^ 



for I„ 



^ 40°). 



APPENDIX B: SLOPE OF THE FLUX 
DISTRIBUTION 

In order to relate the number of resolved point sources to 
the integrated fiux from a population as in the previous 
section, it is necessary to know the slope /3 for the dN/dS 
fiux distribution (eq. (All). In Figure [2] we plot slopes in 
the log A'^ — log 5* panels for the ordinary pulsar and MSP 
models to guide the eye; based on these, we approximated 
the high-fiux end of the MSP distribution to have an 
Euclidean slope (5 — 2.5 in our analytical estimates. In this 
section, we show how the different slopes arise and argue 
that they are robust to the details of our population models. 

Let us first assume that all the sources have an in- 
trinsic luminosity L. We consider three idealized cases, 
corresponding to homogeneous spatial distributions along 
a one-dimension line (linear), in a two-dimensional plane 
(planar), and in three-dimensional space (Euclidean). In 
general. 



dN 

Is 



dN 
dr 



dr 
dS 



(Bl) 



where r is the distance from the observer. Since S — L/iirr'^, 
\dr/dS\ = 2-Kr^ /L for all geometries. If we denote by A, 
(J, and p the linear, planar, and Euclidean spatial source 
densities we have 



dN 
dr 



X 

2nra 
itrr^ p 



if linear 
if planar 
if Euclidean 



Factoring out the constants, this gives 

f f{L)S'''*^^ if linear 



dN 
dS 



9(L)S-^ 
h{L)S- 



/2 



if planar 
if Euclidean 



(B2) 



(B3) 



for /, g, and h that are functions of L and where we have 
written (dN/dS) (L) to emphasize we have assumed that all 
the sources have the same luminosity. 

Consider now the case of an arbitrary luminosity function 
0(1/), normalized such that J dL(j>{L) — 1. Then the total 
fiux distribution is simply given by 

-=/rfL<^(L)-(L). (B4) 



Explicitly, 
dN 

Is ' 



{J dLf{L))S~''^/'^ if linear 
{JdLg{L))S-^ if planar , (B5) 
{J dLh{L))S~^''^ if Euclidean 



which shows that the slope of the distribution is preserved 
regardless of the luminosity function. The dN/ dS slope how- 
ever depends on the geometrical distribution of the sources. 
In cases in which the spatial source distribution is bounded, 
for example, one can show that the luminosity function does 
shape the fiux distribution to some extent. Empirically, we 
find that the high-fiux end of the MSP log A' - log S distri- 
bution has a nearly Euclidean slope (5 = 2.5 (Figure [2| re- 
gardless of the details of the population parameters (which 
do affect the luminosity function). This refiects the near- 
Euclidean nature of the spatial distribution in the Solar 
neighborhood. The distribution departs from this slope at 
lower fiuxes principally owing to the finite scale height of the 
distribution above the Galactic plane. On the other hand, 
the high-flux end of the ordinary pulsar population is inter- 
mediate between the planar (/3 = 2) and linear (/3 = 3/2) 
cases, reflecting the fact that the brightest 7~ray pulsars 
are the young ones located in the Galactic plane and con- 
centrated in spiral arms. 
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164 


0.693 


MSPl_Pbreakl5 


5 


1 


(1.0,1.5,0,-1) 


(8.5, 0.2) 


(0.1,*,-1,*) 


0.0295 


823 


3.76 


MSPl_Pbreak3 


5 


1 


(1.0,3,0,-1) 


(8.5, 0.2) 


(0.1,*,-1,*) 


0.0259 


500 


2.08 


MSPl_Pbreak5 


5 


1 


(1.0,5,0,-1) 


(8.5, 0.2) 


(0.1,*,-!,*) 


0.0228 


342 


1.51 


MSPl_PbreaklO 


5 


1 


(1.0,10,0,-1) 


(8.5, 0.2) 


(0.1,*,-1,*) 


0.0182 


203 


0.938 


MSPl_logB8 


5 


1 


(1.0,*,-1,*) 


(8, 0.2) 


(0.1,*,-1,*) 


0.0347 


84 


0.602 


MSPl_logB825 


5 


1 


(1.0,*,-1,*) 


(8.25, 0.2) 


(0.1,*,-1,*) 


0.0326 


367 


1.54 


MSPl_logB875 


5 


1 


(1.0,*,-!,*) 


(8.75, 0.2) 


(0.1,*,-!,*) 


0.0313 


2320 


11.70 


MSPl_sigB03 


5 


1 


(1.0,*,-1,*) 


(8.5, 0.3) 


(0.1,*,-1,*) 


0.0317 


1275 


6.21 


MSPl_sigB04 


5 


1 


(1.0,*,-1,*) 


(8.5, 0.4) 


(0.1,*,-1,*) 


0.0325 


1605 


9.31 


MSPl_sigB05 


5 


1 


(1.0,*,-1,*) 


(8.5, 0.5) 


(0.1,*,-!,*) 


0.0314 


2148 


12.4 


MSP2_base 


5 


1 


(1.5,*,-1,*) 


(8, 0.2) 


(0.1,*,-1,*) 


0.0245 


4 


0.0775 


MSP2_Lrlow0005 


5 


1 


(1.5,*,-1,*) 


(8, 0.2) 


(0.005,*,-l,*) 


0.495 


195 


1.51 


MSP2_Lrlow001 


5 


1 


(1.5,*,-1,*) 


(8, 0.2) 


(0.01,*,-1,*) 


0.214 


79 


0.698 


MSP2_Lrlow002 


5 


1 


(1.5,*,-!,*) 


(8, 0.2) 


(0.02,*,-!,*) 


0.115 


44 


0.438 


MSP2_Lrlow005 


5 


1 


(1.5,*,-!,*) 


(8, 0.2) 


(0.05,*,-!,*) 


0.0459 


13 


0.175 



Lengths are in kpc, periods in s, magnetic fields in G, radio luminosities in mjy kpc^ at 1.4 GHz, and the 7— ray background 
intensity is in 10~^ ph s~^ cm~^ sr~^. Asterisks indicate that a single power-law index is assumed throughout, with no break. 



